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Abstract 

We discuss the radio emission from high-energy cosmic-ray induced air showers hitting Earth’s surface before the cascade 
has died out in the atmosphere. The induced emission gives rise to a radio signal which should be detectable in the 
currently operating Askaryan radio detectors built to search for the GZK neutrino flux in ice. The in-air emission, the 
in-ice emission, as well as a new component, the coherent transition radiation when the particle bunch crosses the air-ice 
boundary, are included in the calculations. 
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1. Introduction 

We calculate the radio emission from cosmic-ray-induced 
air showers as a possible (background) signal for the Askaryan 
radio-detection experiments currently operating at Antarc¬ 
tica [D HI [3]. A high-energy neutrino interacting in a 
medium like (moon)-rock, ice, or air will induce a high- 
energy particle cascade. In 1962 Askaryan predicted that 
during the development of such a cascade a net negative 
charge excess arises mainly due to Compton scattering [4] . 
This net excess charge by itself will induce a radio signal 
that can be used to measure the original neutrino. This 
Askaryan radio emission BUM has been confirmed ex¬ 
perimentally at SLAG [3 , and more recently the Askaryan 
effect was also confirmed in nature by the radio emission 
from cosmic-ray induced air showers [a Ellin!. 

For high-energy cosmic-ray air showers, along with the 
Askaryan emission, there is another emission mechanism 
due to a net transverse current that is induced in the 
shower front by Earth’s magnetic field [11-14]. Recently 
the radio emission from cosmic-ray air showers has been 
measured in great detail by the LOFAR collaboration nni 
mms], confirming the predictions from several indepen¬ 
dent radio emission models [17-20]. 

Most Askaryan radio detectors [1-3,21-23] search for 
so-called GZK neutrinos that are expected from the in¬ 
teraction of ultra-high-energy cosmic-ray protons with the 
cosmic microwave background The expected GZK 

neutrinos are extremely energetic with energies in the EeV 
range, while the flux at these energies is expected to fall be¬ 
low one neutrino interaction per cubic kilometer of ice per 
year. Therefore, to detect these neutrinos an extremely 
large detection volume, even larger than the cubic kilo¬ 
meter currently covered by the IceCube experiment, is 
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needed. Due to its long attenuation length, the induced 
radio signal is an excellent means to detect these GZK 
neutrinos. This has led to the development of several 
radio-detection experiments [1-6,26-30]. Nevertheless, the 
highest-energy neutrinos detected so-far are those observed 
recently by the IceCube collaboration [3T] and have ener¬ 
gies up to several PeV, just below the energies expected 
from the GZK neutrino flux. 

In this article we calculate the radio emission from 
cosmic-ray-induced air showers as a possible (background) 
signal for the Askaryan radio-detection experiments cur¬ 
rently operating at Antarctica mw- Besides the emis¬ 
sion during the cascade development also transition radi¬ 
ation should be expected when the cosmic ray air shower 
hits Earth’s surface [3H [33] . It follows that the induced 
emission is very hard to distinguish from the direct Askaryan 
emission from a high-energy neutrino induced cascade in 
a dense medium such as ice. 


2. Radio emission from a particle cascade 


We start from the Lienard-Wiechert potentials for a 
point-like four current from classical electrodynamics and 
closely follow the macroscopic MGMR [33] and EVA [5D| 
models. Both models were developed to describe the ra¬ 
dio emission from cosmic-ray-induced air showers. The 
Lienard-Wiechert potentials for a point charge, Ap^{t,x), 
as seen by an observer positioned at x at an observer time t 
are obtained directly from Maxwell’s equations after fixing 
the Lorenz gauge [35] . 
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The point-like current is defined by = eV^, where e 
is the charge, and is the four-velocity for a particle at 
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where the retarded emission time is denoted by tr- 
The denominator of the vector potential, T>, is the retarded 
four-distance. For an extended current with longitudinal 
dimension h and lateral dimensions r, the vector potential 
has to be convolved with the charge distribution given by 
the weight function w{f, h), 
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where the vector potential has to be evaluated at the re¬ 
tarded emission time tj.. The corresponding geometry is 
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Figure 1: The geometry used to calculate the radiation 
emitted from a charge cloud crossing a boundary at z = 
Zb- The observer is positioned at an impact parameter 
d = ^/ix- ryf A{y- TyY. 

denoted in Fig. We consider an observer positioned at 
an impact parameter d = \/{x — A- {y — VyY perpen¬ 
dicular to the charge track, where r^, and Xy denote the 
lateral position of the considered charge within the charge 
cloud. Defining the element in the plane of the observer 
perpendicular to the charge-track as z = 0, we can define 
the time at which the front of the charge cloud crosses this 
plane to be t = 0. Using these definitions the position of 
the charge along the track is now given by z = —ctr + h. 

Fixing the geometry, the vector potential can now be 
evaluated. The retarded emission time is obtained from 
the light-cone condition with respect to the optical path 
length L, 

c{t — tr) = L, (3) 


where the distance di , the distance covered by the emission 
in layer i, is obtained by using a ray-tracing procedure 
based on Snell’s law. Following [35], the retarded distance 
for a signal traveling through different media is given by. 


„ ^ dt 

V = L— . 
dtr 


(5) 


In this work the index of refraction is assumed to be in¬ 
dependent of frequency within the radio frequency range 
starting from a few MHz, up to several GHz. In the simpli¬ 
fied situation where the signal travels through a medium 
with constant index of refraction n, the retarded distance 
can be written in the more common form, 


V = nR{l — n[i cos{9)) , (6) 

where 9 denotes the opening angle between the line of sight 
from the emission point to the observer and the direction 
of movement of the emitting charge. 


2.1. Cherenkov effects for a single electron 

For a single electron moving at a highly relativistic 
velocity fi = vjc k. \ along the z-axis (by definition), 
the current is given by = e (1,0, 0, —/?). The electric 
field is now obtained directly from the Lienard-Wiechert 
potentials through. 
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where i = x,y gives the polarization of the field in the 
transverse direction, and a;® denotes the observer position 
in the transverse plane {x^ = x,x^ = y). For the moment 
we will ignore the electric field in the longitudinal direction 
and, since oc J® = 0 for i = 1, 2 (there is no transverse 
current), we only have to consider the spatial derivative of 
the scalar potential. The electric field in the longitudinal 
direction will in general be small and can easily be calcu¬ 
lated following the gauge condition k ■ e = 0, where k is 
the momentum vector of the photon and e the polarization. 
Hence the photon cannot be polarized along its direction 
of motion. Starting at the zeroth component of the vector 
potential, the spatial derivative can be evaluated by. 
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which corresponds to the radiation from a net charge mov¬ 
ing through the medium. For a relativistic electron (9 ~ 1) 
moving in a medium with a refractive index n > 1 this 
term becomes. 


from which the relation between the observer time and the 
emission time, tr(t), can be obtained. It should be noted 
that tr is a negative quantity. For a medium consisting 
out of m layers with different index of refraction Ui, the 
optical path length can be defined by 

m 

L^'^mdi, (4) 
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Where the label ’st’, denotes that the field is due to a 
highly relativistic non time-varying steady charge. The 
emission shows a radial polarization direction and vanishes 
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linearly with the distance of the observer to the shower 
core. This component of the electric field is suppressed by 
the factor 1 — which vanishes in vacuum. In a medium 
with an index of refraction larger than unity however, this 
factor does not vanish and Cherenkov radiation is observed 
at the point where the retarded distance vanishes, V = 

+ (1 — + j/2) = 0. 

The retarded distance vanishes at the finite Cherenkov 
angle cos{9ch) = (see Eq. ( j ^) where the electric field 
diverges. One intuitive way to understand the Cherenkov 
effect follows from the more general definition of V given 
in Eq. (§. For a vanishing retarded distance, the deriva¬ 
tive dt/dtr has to vanish. It follows that the function t{tr) 
is flat at this point. Hence at an observer time t, signals 
emitted at different emission times will be observed at 
once, leading to a boosted electric field. The vanishing 
of the retarded distance leads to a divergence in the elec¬ 
tric field expressions. These divergences are integrable and 
therefore disappear for coherent emission by performing an 
integration over the finite charge and current distributions 
in the shower front [20] . 


2.2. Transition radiation for a single electron 

So far we calculated the component of the electric field 
due to a relativistically moving net charge in a medium 
with a refractive index equal to n. How does this compare 
to the transition radiation for a relativistic charge cross¬ 
ing from a medium with refractive index ni to a medium 
with refractive index n 2 ? The vector potential for a single 
electron now becomes. 


where the discontinuity at a distance zi, = —cti,, corre¬ 
sponding to the retarded emission time tb when the elec¬ 
tron crosses the boundary, is reflected by the step function 
9{x) which is defined by, 



if a; < 0 
if a; > 0 


( 11 ) 


Since the step function is a function of the retarded emis¬ 
sion time. 


9{z - Zb) = e{-c{tr - tb)) , (12) 

an additional term has to be added to Eq. The full 
electric field is now given by. 
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where in case of a single electron the second term on the 
right hand side of Eq. (131 will correspond to the transi¬ 
tion radiation. The transition radiation can therefore be 


evaluated as. 
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It follows that when there is no boundary, hence ni = n 2 , 
the transition radiation vanishes as it should. Looking 
more closely at the obtained expression in Eq. (14), it can 
be described as the superposition of the emission just be¬ 
fore the particle crosses the boundary and the field just 
after the particle crossed the boundary. The two terms in¬ 
terfere destructively. This corresponds well to the mirror- 
charge approach for determining the transition radiation 
as applied by Ginzburg et al. m and the expressions ob¬ 
tained in ESI ISSl and references therein. 



Figure 2: The angular intensity of the transition radiation 
seen from a charge moving from vacuum, ni = 1, to a 
dense medium n 2 > 1. 


In Fig. the intensity of the transition radiation is 
shown as function of angle. Since in vacuum the retarded 
distance vanishes at forward angles the first term domi¬ 
nates in the forward direction. At very small angles the 
intensity is suppressed due to the factor x in Eq. (14). 
With increasing angle l/\'D\f drops in magnitude while 
the contribution from below the surface, proportional to 
l/lT’lt _e increases in magnitude. Since the two interfere 
destructively there is a cancellation at a certain angle af¬ 
ter which the second term dominates. At the Cherenkov 
angle in the medium, 02 j the second term diverges. 


2.2.1. The sudden appearance signal 

Another effect which is similar to transition radiation 
is the sudden appearance signal. One example of such a 
signal is found in accelerator experiments [39l |40]. Here 
the charged particle beam is accelerated, but during the 
acceleration process the charge is (partly) blocked from an 
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observer. The consequence is that when the beam leaves 
the accelerator, the observer suddenly observes a charge. 
This effect can be described in a similar way as transition 
radiation. The vector potential simply becomes, 

A°{t,x)= -^-^9{zb-z), (15) 


for a beam traveling in a medium with refractive index n. 
Following the transition radiation calculation, the electric 
field is now directly obtained by, 
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It should be noted that the delta-function is a function of 
the retarded emission time, tr- The functional dependence 
can be shifted to the observer time, t, after which the field 
is given by the more common expression, 


the radio emission arising from the transverse current the 
reader is referred to [501131] • 

In case of a point-like time varying charge Q{tr) = 
eNf.{tr), the partial derivative with respect to the retarded 
time of the second term on the right hand side of Eq. (131 
will now also get a contribution where the derivative acts 
on this current. This leads to the varying charge emission, 
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which is the far-field radiation component. It should be 
noted that also here the signal scales like I/IT’P and hence 
Cherenkov or equivalently relativistic time-compression ef¬ 
fects apply equally well for this component of the radia¬ 
tion. Furthermore, a similar polarization behavior as for 
the emission from a highly relativistic steady charge as 
well as the transition radiation is obtained. 


3. Coherent emission 
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where Lb denotes the optical path length for the signal 
emitted at the boundary point toward the observer. 


2.3. Time varying current emission 

So far we considered radiation from a single electron. 
In case of an electron bunch, there will be another radi¬ 
ation component due to the time variation of the total 
number of charges. In general this time variation can be 
linked to the net contribution of coherent bremsstrahlung 
emission of charges dropping out of the high-energy charge 
cloud and the emission of Compton electrons which are 
suddenly accelerated to relativistic speed. More generally, 
we can define the total number of particles at the retarded 
emission time tr by the distribution N^itr). Defining the 
four-current as. 


J^^{tr) = eNe{tr)V^^ , 


the vector potential becomes. 
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which gives the point-like vector potential for a non ex¬ 
tended current. For a cosmic-ray air shower, the two main 
emission mechanisms are due to a time-varying transverse 
current which is induced by Earth’s magnetic field, and 
the Askaryan emission due to the time-variation of the 
net negative charge-excess in the cascade. Nevertheless, 
for the geometry considered in this article, describing the 
emission for a perpendicular incoming shower hitting the 
ice surface at the South-Pole, the cascade will be aligned 
with Earth’s magnetic field and the transverse current van¬ 
ishes. Therefore, in this section we focus on the emission 
from a time varying charge. For more information about 


To obtain the coherent emission we need to consider 
the spatial extent of the particle cascade. This is done 
by inclusion of the weight function w{f,h). The weight 
function is normalized such that f dh d^rw(r,h) = 1. In 
a realistic situation, there will also be emission from the 
charged trail which is left behind after the cascade has 
passed. A detailed calculation including this can be found 
in |34j . The expressions given below for the coherent emis¬ 
sion are obtained including this positive trail. 

The electric field is now obtained by convolving the 
point-like current with the particle distributions in the 
shower front which can be evaluated numerically, 
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Here p = eV x (eV x C/s) is the polarization of the signal, 
is the unit vector pointing from the emission point to the 
observer, and ep is the unit vector denoting the direction of 
the cascade. For the transition radiation the delta-function 
in Eq. ( |l4| ) can be rewritten as, 

5{z — Zb) = S{h — c{tr — tb)) . (23) 


The electric field can now be solved analytically by inte¬ 
grating the delta-function and is given by. 
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h = c{tr — tb) 


Following the same procedure the sudden appearance sig¬ 
nal is given by, 


In the previous section, we obtained the electric field 
expressions for the transition radiation from a particle cas¬ 
cade traversing the boundary between two different media. 
We also considered the steady charge emission as well as 
the varying current emission. We now have all ingredients 
to solve for the the emission from a high-energy cosmic-ray 
air shower which penetrates a surface. 
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(25) 4-1- The particle cascade 


h — c{tr — tb) 


4. The cosmic-ray air shower signal in Askaryan 
radio detectors 
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Figure 3: The number of charges W as a function of 
height z above the antenna as simulated for a cosmic-ray 
air shower with primary particle energy of 10^^ eV. The 
antenna is positioned at an observer level 2900 m above 
sea level in ice. The ice-air boundary is at 3000 m above 
sea level. 


As an example we consider a shower which is induced 
by a 10^^ eV primary cosmic ray, where the shower will hit 
an ice surface. The shower profile can be expressed as a 
function of the penetration depth X{g/cm?') = f p(z)dz, 
given by the line integral over the density which the shower 
has passed through. This allows us to naturally take into 
account for the air-ice boundary by simply writing the 
density as. 


p(z) = Pair 0{z - Zb) + Pice 0{zb - z) , 


(26) 


again using the step function 9{x). We will assume a 
density profile Pair{z) = poexp[— (^(z + zq)] for an ob¬ 
server positioned at a height z = zq above sea level. Here 
C = 1.160 X 10“"^ m“^, and po = 1.168 x 10“^ g cm“^ [M]. 
The ice density is assumed to be constant over the few 
meters in which the cascade will die out and taken as 
Pice = 0.92 g cm“3. 

Since the radiation length Xq = 36 g/cm^, as well as 
the critical energy Ecrit = 80 MeV for electrons is ap¬ 
proximately equal in air and ice, we can now take a NKG 
approximation SDIU] given by. 


N{X) = 


0.31 exp[(A/Ao)(l - 1.51ns)] 
iy\n{E / Ecrit) 


(27) 


for the total number of particles as a function of depth. 
The shower age s is given by [33], 


s(A) = 


3X/Xo 


{X/Xo) + 2HE/Ecrit) 


(28) 


The excess charge as function of shower depth can be 
approximated by Nch{X) « 0.23 A(A) [34]. The total 
number of excess electrons as function of depth is shown 
in Fig. Taking a geometry with the observer positioned 
at 2900 m above sea level with the air-ice boundary at 
3000 m, the boundary at Zf, = 100 m is clearly visible. 

The particle distribution in the shower front is given by 
the weight function u>(r, h) = 6{r)f{h), where / dh f{h) = 
1. The radial extension of the particles in the shower front 
is taken to be a delta function at the shower axis. To com¬ 
pensate for the loss of the lateral coherence scale, we use an 
effective width hi for the longitudinal particle distribution 
in the shower front. This width therefore reflects the co¬ 
herence scale due to the full extension of the shower front. 
The longitudinal distribution of particles in the shower 
front is parametrized by [331 |33j , 

f{h) = {A/hl)he-^'^/'^^ . 
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The effective width hi is chosen to be hi" = 0.5 m fol¬ 
lowing [45] for the in-air development, while for the in-ice 
part of the cascade a width h"^ = 0.1 m is chosen. 

4-2. The refractive index 

In [45j it was shown that in determining the radio sig¬ 
nal from cosmic-ray air showers it is crucial to take into 
account a realistic index of refraction. Therefore, in the 
following we model the index of refraction in air by the 
Gladstone-Dale law, 

nair{z) = 1 -1- 0.226-^Pairiz) . (30) 

cm'^ 

Furthermore, in [36j it was shown that the bending of the 
emission in air can safely be neglected. The index of re¬ 
fraction in ice is taken as a constant equal to 

nice = 1.78 . (31) 


4-3. Results 

Since we now have our electric field expressions, as well 
as the particle distributions, we can calculate the electric 
field at different observer positions in ice. We consider two 
different boundary levels at 500 m and 3000 m above sea 
level. The shower profile is given in Fig. for a geome¬ 
try where the air-ice boundary is 3000 m above sea level. 
In Fig. 3a the full shower profile is given. It follows that 
the shower hits the ice surface before it reaches its maxi¬ 
mum. This is more clear from Fig. |3b| where we zoom in 
on the boundary. In the ice the shower quickly reaches its 
maximum and dies out within 10 meters. 

We consider the emission as seen by an observer po¬ 
sitioned at several distances, c? = 40 m, d = 80 m, and 
d = 240 m, from the shower axis, 100 m below the air-ice 
boundary. The obtained electric fields are shown in Fig. 

a, c, and e. The full red line gives the electric field gener¬ 
ated by the in-air development of the cascade, the striped 
purple line gives the transition radiation, and the dotted 
blue line gives the emission from the in-ice development of 
the cascade. Next to the obtained electric fields, in Fig. 

b, d, and f we also show the total number of particles 
as function of height. Furthermore, in these figures we 
show the observer time for a signal emitted from a cer¬ 
tain height. The full green line gives the total number of 
particles as function of height. It should be noted that 
the vertical axis is shifted by 90 m for plotting purposes. 
The full red line gives the emission height as function of 
the observer time for the in-air emission, where the striped 
purple line gives the same quantity for the in-ice emission. 

The emission observed at d = 40 m is shown in Fig.|^. 
As follows from Fig. [^, the in-air emission from large 
heights is observed before the emission from lower heights. 
For the in-ice emission, this picture is completely reversed. 
In-ice signals emitted from large heights are delayed by the 
medium, while the cascade continues to propagate with 
the speed of light. Hence signals emitted at later times 
(lower heights) arrive before signals emitted early and the 


observer is positioned inside the Cherenkov cone for the 
full in-ice emission. 

For the transition radiation, it is important to notice 
that the electric field as given in Eq. (14) can be seen 
as a superposition of the emission just above, and just 
below the boundary which interferes destructively. The 
emission scales with 1/|I?| \dtr/dt\ ~ \dz/dt\, which is 
reflected in Fig. & by the slopes of the full red and striped 
purple lines at the boundary. At the boundary, there will 
be a sudden change of the particle distributions in the 
shower front. To take this into account in our modeling, 
the emission just above the boundary is evaluated using 
the particle distributions for the in-air shower, where the 
component just below the boundary is modeled using the 
particle distribution for the in-ice cascade. From this point 
of view one might also consider the transition radiation 
from just above the boundary as the sudden-death signal 
from the emission in air, where the transition radiation 
from just below the boundary can be seen as the sudden 
appearance signal for the in-ice emission. 

For an observer positioned at d = 80 m, see Fig. H: 
and d, a similar picture is obtained. The in-air emission 
is observed over a longer time-scale since we move further 
away from the Cherenkov cone for the emission emitted 
in air. For the in-ice emission, however, we shift closer to 
the Cherenkov angle. It follows that the in-ice emission 
is observed within a much shorter time-span and becomes 
much stronger. The transition radiation is now dominated 
by the emission from just below the boundary. 

Finally we consider an observer positioned at d = 240 m. 
From Fig. and f it follows that both the in-air emis¬ 
sion as well as the in-ice emission are observed outside the 
Cherenkov cone. The emission is observed over a rather 
long time-scale, although the in-air component starts to 
be rather weak. One interesting feature is that the emis¬ 
sion just above the boundary does not arrive at the same 
time as the emission just below the boundary. The tran¬ 
sition radiation component just above the boundary is 
highly suppressed, and arrives just before t = 1000 ns, 
where the component emitted just below the boundary is 
much stronger and arrives at a much later time around 
t = 1200 ns. The time difference arises due to the fact 
that the signal emitted just above the air-ice boundary 
will first travel a significant part of its path almost hor¬ 
izontally before breaking into the ice under the critical 
angle, which in this case is equal to the in-ice Cherenkov 
angle. The signal emitted just below the air-ice boundary 
will travel its full path through the ice and hence obtain 
a large delay with respect to the signal emitted just above 
the boundary. It should be noted that this effect occurs 
in the situation of a perfectly flat and smooth surface. In 
a realistic experiment, the emission from just above the 
boundary however will not be able to travel almost per¬ 
fectly horizontally along the surface and hence will loose 
coherence and become suppressed (already in the present 
case it is almost negligible in magnitude). The signal emit¬ 
ted just below the surface will not be affected and keep its 
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Figure 4: The electric field at different observer distances equal to, a) c? = 40 m, b) d = 80 m, c) d = 240 m. The figures 
on the right show the emission height, plotted as function of the observer time. The full red line gives the emission in air, 
the dotted purple line gives the transition radiation, and the dashed blue line gives the in-ice emission. For the figures 
on the right, the total number of particles is given by the f^l green line. 



































coherence. 

The emission will be coherent up to relatively high 
frequencies. This is also seen in Fig. 5a where we plot 
the frequency spectrum of the different components of the 
emission when the observer is positioned at a distance of 
d = 240 m. In Fig.[^ we plot the frequency spectrum for 
the same geometry shifting the ice-air boundary to 500 m 
above sea level. Coherence of the in-ice emission as well 


Electric field (Z(,=3000 m) 



(a) 


Electric field (zi,=500 m) 



(b) 

Figure 5: The frequency spectrum of the different compo¬ 
nents to the electric field as seen by an observer positioned 
100 m below the ice-air boundary and d = 250 m from the 
shower axis. The simulation is performed for a 10^^ eV 
primary energy cosmic-ray air shower. 

as the in-air emission away from the Cherenkov angle is 
typically determined by the length of the shower trajec¬ 
tory leading to a suppression at the highest frequencies. 
The transition radiation, however, is emitted from a single 
point at the boundary, and hence its coherence is fully de¬ 
termined by the particle distributions in the shower front 


which gives a cut-off at relatively high-frequencies in the 
GHz range. Each of the several different components has 
a finite response at zero frequency. One should note how¬ 
ever that the combined response of all different compo¬ 
nents vanishes at zero frequency. 

In Fig.|^ and Fig.[^ the detailed properties of the emis¬ 
sion in time and frequency space were shown. This allows 
us to understand the angular distribution of the different 
components of the emission shown in Fig. Here we plot 
the integrated absolute value of electric field I = J \E\dt. 
It follows that the in-air emission given by the full red 




(b) 

Figure 6: The angular distribution of the integrated field 
for a radial observer distance of i? = y/z^~+~d^ = 300 m. 
The full red line gives the angular distribution for the in¬ 
air emission, the dotted purple line gives the transition 
radiation, and the dashed blue line shows the in-ice emis¬ 
sion. 

line in Fig. as expected, peaks toward a highly forward 
angle {9 = 0°), after which the emission drops rapidly to- 
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ward larger angles. The transition radiation, shown by 
the dotted purple line in Fig. shows a similar behav¬ 
ior as observed in Fig. There is a strong forward peak, 
after which the transition radiation gets suppressed due 
to the destructive interference between the emission just 
above the air-ice boundary and the emission just below 
the air-ice boundary. When the destructive interference 
is maximal, contrary to Fig. the field does not vanish 
completely. This is due to the different particle distribu¬ 
tions for the in-air emission and the in-ice emission which 
is taken into account for in the evaluation of the transi¬ 
tion radiation. The in-ice emission, as expected, peaks 
at the in-ice Cherenkov angle. Nevertheless, the emission 
pattern is rather broad toward smaller angles due to the 
longitudinal extent of the cascade. 

4.4- Zenith angle dependence 

Up to now we only considered a perpendicular incom¬ 
ing cosmic-ray induced air shower. Since a shower coming 
in under a finite zenith angle can be treated by a direct 
rotation of the geometry, we do not expect the emission 
to change significantly. One effect that is to be expected 
for an inclined shower, or in case of a non-perpendicular 
boundary, is that the transition radiation from different 
radial parts of the shower is emitted at different times. 
In case of an observer positioned underneath the shower 
axis this will lead to a small additional spread in the ar¬ 
rival time of signals emitted from different positions in the 
shower front, and hence a slight decrease of pulse-strength 
can be expected. For an observer positioned away from 
the shower axis however, this effect is reversed, leading to 
a slightly enhanced pulse-strength in the detector. Fur¬ 
thermore, it should also be noted that for more inclined 
geometries, a larger part of the signal created in air will be 
reflected off of the surface suppressing the in-air emission 
over the in-ice emission even more. 

The most important effect, however, will be due to 
the change of the total number of charges hitting the air- 
ice boundary. Since for larger zenith angles the shower 
traverses a longer distance through air, the total number 
of particles actually hitting the air-ice boundary changes. 
Other effects influencing the total number of particles hit¬ 
ting the air-ice surface are the chemical composition and 
the energy of the primary cosmic ray. An iron induced 
shower typically develops earlier in the atmosphere than a 
proton induced shower, where cosmic rays of higher energy 
typically peak deeper in the atmosphere. In Fig. [7] we show 
the total number of particles hitting the air-ice surface for 
boundary layers at Zb = 3000 m (full lines) and Zb = 500 m 
(dotted lines) for a typical proton shower with a primary 
energy of Ep = 10^^ eV (red lines) and Ep = 10^® eV (blue 
lines). 

For an air-ice boundary at Zb = 3000 m, the air shower 
is still below shower maximum for both considered ener¬ 
gies. It follows that the total number of particles peaks at a 
zenith angle of approximately 0 « 40 — 50 degrees, where 
the air shower is fully developed at the boundary. For 


larger zenith angles, the total number of particles hitting 
the air-ice boundary becomes smaller, and for zenith an¬ 
gles larger than 0 > 60 degrees the shower dies out before 
hitting the air-ice boundary. Hence no transition radiation 
and in-ice emission will be observed for showers at zenith 
angles larger than approximately 0 « 60 degrees. 



Figure 7: The total number of particles hitting the air- 
ice surface for boundary layers at Zb = 3000 m (full lines) 
and Zb = 500 m (dotted lines) for a typical proton shower 
with a primary energy of Ep = 10^^ eV (red lines) and 
Ep = 10^® eV (blue lines). 

4-5. Cosmic-ray air shower or neutrino induced cascade? 

One important question to consider is how the cosmic- 
ray air shower signal compares to the emission from a neu¬ 
trino induced cascade in ice. This is shown in Fig. For 
the cosmic-ray air shower signal, we consider both the in¬ 
ice emission as well as the transition radiation component 
just below the boundary. As follows from Fig. the in-air 
emission is very small and will therefore be ignored for this 
comparison. 

The neutrino induced cascade is modeled by taking a 
10^^ eV primary neutrino interacting at the surface of the 
air-ice boundary, after which the cascade develops in ice. 
The observer is positioned 100 m below the ice surface at a 
distance d = 240 m perpendicular to the shower axis. The 
effective width of the particle distribution is taken as hi = 
0.1 m, approximately corresponding to the dimensions of 
the cascade front in ice. 

It follows that the emission from a typical proton in¬ 
duced high-energy cosmic ray air shower hitting the ice 
surface is of similar magnitude compared to the Askaryan 
signal from a neutrino induced particle cascade in ice of 
the same energy. Since the total number of particles hit¬ 
ting the air-ice interface is slightly larger at Zb = 500 m for 
the considered geometry of a vertical cosmic-ray air shower 
(see Fig.[^, the transition radiation gets slightly enhanced 
with respect to the emission seen at Zb = 3000 m. 
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Figure 8: The electric field from a neutrino induced par¬ 
ticle cascade in ice (full red line), compared to the transi¬ 
tion radiation just below the boundary layer and the in¬ 
ice emission from a vertical, proton induced, cosmic-ray 
air shower of the same energy. The electric field is shown 
for ice-air boundaries equal to Zb = 3000 m (dashed blue 
line) and Zf, = 500 m (dotted purple line). The observer is 
positioned at a lateral distance d = 250 m perpendicular 
to the shower axis. 


The direct consequence is that the radio emission from 
a high-energy cosmic ray air shower will be very hard to 
distinguish from the emission of a high-energy neutrino 
indnced particle cascade without directional information 
about the signal. Another possible method to separate the 
cosmic-ray air shower from the neutrino induced cascade 
might be obtained by adding a surface veto to the array. 

5. Summary and conclusions 

We derived an analytical expression for coherent tran¬ 
sition radiation from a particle bunch with a net charge 
which is traversing from one medium to another. In ad¬ 
dition to the transition radiation also the emission due to 
a highly-relativistic steady charge and a time-varying cur¬ 
rent are given. 

As a first application we calculate the radio emission 
from a cosmic-ray-induced air shower hitting an ice layer 
before the shower has died out. It is shown that a rel¬ 
atively strong transition radiation component can be ex¬ 
pected from an air shower when taking an air-ice bound¬ 
ary in the range between 500-3000 m above sea level. The 
emission from such a cosmic-ray air shower is calculated 
to be of similar strength as the Askaryan signal obtained 
from an in-ice cascade induced by a GZK-neutrino of sim¬ 
ilar energy. Furthermore, the polarization of the transi¬ 
tion radiation will be similar to the polarization of the 
Askaryan signal. It follows that without directional infor¬ 
mation or a surface veto, it will be very hard to distin¬ 
guish between both signals. Therefore, the emission from 


cosmic-ray air showers might induce a significant (back¬ 
ground) signal in the currently operating Askaryan neu¬ 
trino detectors at Antarctica. 

Next to the discussed application for cosmic-ray air 
showers hitting a dense medium, the transition radiation 
from neutrino induced particle cascades traveling from a 
dense medium to air or vacuum is also expected to give 
a strong signal. This signal might be a promising probe 
to detect high-energy neutrino-induced particle cascades 
escaping dense media. A more detailed calculation for this 
component will be given in future work. 
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